Stability of RVB hole stripes in high-temperature superconductors 
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Indications of density-wave states in underdoped cuprates, coming from recent STM (scanning 
tunneling microscopy) and Hall-resistance measurements, have raised new concerns whether stripes 
could be stabilized in the superconducting phase of cuprate materials, even in the absence of an- 
tiferromagnetism. Here, we investigate this issue using state-of-the-art quantum Monte Carlo cal- 
culations of a t — J model. In particular we consider the stability of unidirectional hole domains 
in a modulated superconducting background, by taking into account the effect of tetragonal-lattice 
distortions, next-nearest neighbor hopping and long-range Coulomb repulsion. 



I. INTRODUCTION 

The physics of high-temperature superconductors 
(HTSC) is characterized by the presence of several com- 
peting orders. Besides antiferromagnetism and super- 
conductivity, it is now widely accepted that the un- 
derdoped region of the HTSC phase diagram shows 
a pseudogap phasei"^ and occasionally the presence of 
spatially-ordered states^ii at special dopings, character- 
ized by charge domain walls separated by antiferromag- 
netic regions. Recently, accurate STM experiments^ 
have revealed modulations of the local density of states 
in the form of unidirectional domains for two different 
underdoped HTSC materials, Ca2-5Na5Cu02Cl2 and 
Bi2Sr2CaCu20g+5, where spatial ordering was not de- 
tected before. These findings highlight the importance of 
inhomogeneous states as a general feature of HTSC, and 
extend the idea of quasi-one dimensional order beyond 
the original spin-stripe picture^, with superconductivity 
coexisting with charge modulation, in the absence of anti- 
ferromagnetism. Moreover, recent measurements of Hall 
resistancoi on underdoped YBa2Cu305 and YBa2Cu408 
give evidence of a reconstruction of the Fermi surface 
caused by the onset of a density- wave phase in the large- 
field induced normal state, indicating that spatial sym- 
metry breaking is a common characteristic of underdoped 
cuprates. 

Theoretically, it is challenging to characterize a mod- 
ulated state that could mimic the experimental findings, 
and to identify the key microscopic parameters that could 
induce this kind of non-uniform superconducting phases. 
Indeed, besides the strong, local repulsion among elec- 
trons, that leads to antiferromagnetism, the low-energy 
physics of cuprates is often determined by additional 
hopping integrals among next-nearest neighbor copper 
atoms, and by a long-range Coulomb repulsion among 
electrons. These terms have been studied in the past for 
different inhomogeneous states and were found to stabi- 
lize modulated superstructures^!^. Also of relevance, is 
the role played by structural instabilities, a feature un- 
doubtedly present in many HTSC compounds, in stabiliz- 
ing modulated superstructures. Indeed, HTSC are often 
characterized by a low-temperature tetragonal phase^*'. 



produced by a tilt of the oxygen octahedra that leads to a 
different electronic hopping and exchange integral along 
the two planar directions. Previous calculations found 
that spatial anisotropy could stabilize antiferromagnetic- 
stripe phaseaiiii^. Therefore, in view of what is found 
in experiments, it is important to understand if a similar 
effect occurs for superconducting modulated states. 



In this paper, we study the competition among ho- 
mogeneous and non-uniform superconducting states, by 
taking into account the presence of spatial lattice distor- 
tions, next-nearest neighbors hopping terms and long- 
range Coulomb repulsion. As a prototype model for 
cuprate superconductors, we investigate the t—J modeli^ 
in two dimensions, by using state-of-the-art Variational 
Monte Carlo techniques (supplemented by mean-field ap- 
proaches). Following Anderson idea^^, the variational 
approach, based on the resonating-valence-bond (RVB) 
wavefunction, has successfully described most of the 
features of HTSC^^'^^. However, to tackle spatially- 
inhomogeneous superconducting states like those found 
by experiments, more involved calculations, including the 
possibility of local modulations in the superconducting 
state, are required. Here, we accomplish this task and 
show that, generically, half-filled hole stripes can form 
easily in a RVB superconductor at quite small energy 
cost. Interestingly enough, lattice distortion and long- 
range Coulomb repulsion are proved to play a key role. 



This paper is organized as follows. In Section |TT] we in- 
troduce the model and the variational approach, present- 
ing the two families of modulated wavefunctions that we 
study, the 7r-phase shift and inphase domain RVB wave- 
functions. In Section Hm we characterize in detail the new 
inphase domain RVB state, introduced in this paper, for 
the t — J model. In Section llVi IVl and IVII we show the 
effect of lattice distortion, next-nearest neighbor hopping 
and long-range Coulomb repulsion, respectively, on the 
stabilization of modulated superconducting states. Fi- 
nally, in Section [VIII we draw our conclusions. 



II. MODEL AND WAVEFUNCTIONS 

The t — J Hamiltonian is defined, using standard no- 
tations, asH; 

Ht-j = XI o^vicL^jcr + h.c.) + J XI ■ 

where = (1 — n,;_(j)c|^ acts on the reduced Hilbert 
space with no double occupancies. For non-distorted lat- 
tices we take aij = ax = oty = 1 for all nearest neighbors. 
Instead, tetragonal distortion is set by taking different 
hopping and exchange parameters along the x and y di- 
rection, according to the position of the tilt axis. This 
implies ax 7^ ay, with < 1 {ay < 1) and Oy — 1 
{ax = 1) respectively for the tilt axis along the y {x) 
direction. Here, we mainly focus on the physically rele- 
vant case i/J = 3 (if not specified otherwise) at doping 
S = 1/8. 

In the following, we make use of the Variational Monte 
Carlo techniques (VMC), using clusters of iV 8 x 8, 
128 (45 degree 8^ + 8^ tihed lattice) and 16 x 16 sites 
with periodic-boundary conditions. Moreover, to have 
access to larger sizes that can approach the thermody- 
namic limit, we compare our variational results with pre- 
dictions coming from the renormalized mean-field theory 
(RMFT)il, on a 128 x 128 cluster, using unit ceU transla- 
tion symmetrjsi^ (for more details on the RMFT results, 
see Reffi^). In the variational procedure, the RVB state 
is derived from the BCS mean-field state by incorporating 
the effect of correlation via the so-called Gutzwiller pro- 
jector Vg, that can be treated exactly within the quan- 
tum Monte Carlo scheme. Indeed, given the mean-field 
Hamiltonian: 

Hmf = ^{XiAaCjcj + h.c) 
+ X^^^J^It^i + '^•^■) + mX"- (2) 

(ij) iff 

we construct the variational state by applying the 
Gutzwiller projector to the ground state \D) of Eq. 

\^G)^Vy\D) = ^{\-n,^n,^)\D) (3) 

i 

where riia counts the number of electrons of spin a on site 
i. The terms Xih -'^l- © correspond to a 

set of parameters that are optimized in order to minimize 
the variational energy, according to the stochastic mini- 
mization algorithni^^ . In the simple case of the uniform 
RVB state, the independent parameters reduce to /i and 
Ay- = ±A for nearest neighbors, with the sign follow- 
ing the d-wave symmetry. Moreover, another competing 
homogeneous state, the staggered- fiux phase^i (SFP), is 
obtained by allowing the hopping parameters Xij to be 
complex, leading to staggered currents circulating in op- 
posite directions in the neighboring plaquettes. 
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FIG. 1: (Color online) Cartoon of the (a) bond-centered 
and (b) site-centered domain RVB (DRVB) stripes. Non- 
equivalent bonds have diflferent colors. Thicker (larger) bonds 
(circles) have larger pairing (hole) densities. Dashed lines in- 
dicate bonds where A^w is set to zero. 



Here, we study possible instabilities of the RVB state 
towards spatial modulations, allowing the RVB bonds 
Aij to become inhomogeneous. In particular, it is known 
from experiments that one-dimensional hole-rich regions 
are generally spaced by 1/(2(5), implying half- filled hole 
domains with an average of one hole every 1/(5 sites in 
the direction where translational symmetry is broken. At 
a doping of (5 = 1/8, this corresponds to a periodicity 
of four lattice spacings 4a. Remarkably, we found such 
states, where hole-rich stripes coexist with modulated su- 
perconducting domains, whose energy per site is only of 
a very small fraction of J higher than its uniform coun- 
terpart. The hierarchy of phases is therefore expected to 
be very sensitive to the details of the microscopic Hamil- 
tonian. 

In a previous paper—, we considered a variational 
domain- RVB state that contains a 7r-phase shift in the su- 
perconducting order parameter between regions of four- 
lattice spacings width (ttDRVB). The phase shifts cre- 
ate domain walls with vanishing pairing amplitude and 
a consequent concentration of holes. The ttDRVB states 
turned out to be competitive with respect to the homo- 
geneous RVB state, and could explain the copper-oxygen 
layer decoupling found in La2-£cBaa;Cu04, accompanied 
by a depression of Tc at 1/8 doping^^'^'*. Here, we com- 
plete the characterization of ttDRVB states, by consider- 
ing the effect of lattice distortion, nearest-neighbor hop- 
ping term, and long-range Coulomb repulsion. More- 
over, since the presence of the antiphase domains implies 
a certain amount of energy, in the following we investi- 
gate an alternative strategy to tackle the problem, and 
propose another possible candidate for superconducting 
stripes, where domains are inphase: the domain RVB 
state (DRVB). Starting from the RVB state, we inves- 
tigate the energy cost associated to the introduction of 
line defects in the uniform state, in the form of unidirec- 
tional hole domains. This is done by imposing a vanish- 
ing pairing amplitude Ay — along one direction, with 
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TABLE L VMC and RMFT energy per site (in units of t) 
for different projected wavefunctions for the t — J model at 
doping 1/8. 



0.14 



[t] Ee 



[t] 



RVB 
SFP 
ttDRVB [Ref.2 
BC-DRVB 
SC-DRVB 



-0.45564(3) 
-0.44630(3) 
-0.44529(3) 
-0.45490(3) 
-0.45525(3) 



-0.4549(1) 
-0.4286(1) 
-0.4413(1) 
-0.4511(1) 
-0.4507(1) 



a periodicity that, at doping 1/8, corresponds to 4a. In 
Figure [T] we show the two possible modulations of the 
inphase DRVB state derived from the above arguments. 
In Figure [H^a) stripes are bond-centered (BC), with the 
RVB bonds characterized by a periodicity of 4a along 
the X direction. The pairing field A^- associated to one 
type of horizontal bonds is set to zero. Alternatively, in 
Figure [Ijb) the stripes are site-centered (SC), where the 
lines of vertical bonds with Ay = are separated by 4a 
along the x direction. 

Besides the 7r-shift domain RVB stripes, which could 
give a natural explanation to the decoupling among 
superconducting-layers, this second class of supercon- 
ducting hole stripes, the DRVB states, could be related 
to the STM experimental findings of inhomogeneous su- 
perconducting states. 
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III. DRVB STRIPES IN THE t - J MODEL 
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In this section we characterize the new inphase domain 
RVB states for the t — J model, showing that, among 
the superconducting stripe state o^^i^^ , they constitute a 
valid candidate for a strongly-competing, inhomogeneous 
state. We start by considering the variational energies of 
the different stripe states, together with the correspond- 
ing RMFT values in Table [H Both approaches show that, 
with respect to the SFP and the ttDRVB^ states, the 
DRVB states are one order of magnitude closer in en- 
ergy to the homogeneous RVB phase. Note that the site- 
centered variational energy approaches the RVB state 
very closely. However, in both cases the small energy 
difference indicates both DRVB states as very promising 
candidates for a competing inhomogeneous state. Re- 
markably, within the DRVB geometry, the suppression 
of some of the RVB bonds does not imply significant 
energy costs. Notice that, at the variational level, the 
optimal wavefunctions preserve the d-wave symmetry of 
the pairing field in the regions where A^ 7^ 0, with a 
small modulation that depends upon the type of bond 
(see Figure [1]). 

To prove that the non-homogeneous DRVB wavefunc- 
tions that we have constructed indeed correspond to spa- 
tially modulated states, we show in Figure [5] typical pro- 
files in the bond- and site-centered states. It turns out 
that the hole distribution along the different sites of the 



FIG. 2: (Color online) Top panel: Hole distribution in DRVB 
states for different values of t/ J as shown in graph. Middle 
and bottom panel: Magnetic (in units of J) and kinetic (in 
units of t) bond energy within the unitary cell of Figure [T] for 
t/J = 3. Squares indicate energies associated to horizontal 
bonds (located between two sites), triangles to vertical bonds. 
Left panels refer to the BC-DRVB state and right panels to 
the SC-DRVB state. All quantities are calculated for a cluster 
of 16x16 sites. 



unitary cell is non-uniform, the hole density being larger 
along the regions where A^ = (notice that the magni- 
tude of the hole modulation increases with J/t). More- 
over, the magnetic and kinetic energies follow the same 
modulation than the charge. Indeed, along the Ay = 
bonds the kinetic energy is enhanced, while the magnetic 
energy is suppressed, as expected. Note that the bonds 
presenting the largest modulation are the horizontal (ver- 
tical) ones for the BC (SC) wavefunction, i.e., the bonds 
along the direction where some of the pairing fields are 
suppressed. 

Finally, to characterize the superconducting properties 
of our variational states, we calculate the singlet pairing 
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FIG. 3: (Color online) Top panel: Pairing order parameter 
(filled points for VMC, empty points for RMFT) for the BC- 
DRVB (left) and SC-DRVB (right) states along the 8 inde- 
pendent sites of Figure [1] Squares correspond to horizontal 
bonds (located between two sites), triangles to vertical bonds. 
The blue line corresponds to the pairing order parameter in 
the homogeneous RVB state. Inset: size scaling of the pairing 
order parameter as a function of the inverse- squared number 
of sites. Colors refer to the corresponding colored bonds of 
Figure [T] Stars correspond to the projected Fermi sea, for 
comparison. Bottom panel: corresponding pairing variational 
parameters for the different bonds. 
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FIG. 4: (Color online) Energy difference between the BC- 
DRVB (SC-DRVB) and the homogeneous RVB stated as 
a function of the lattice distortion Ux/oy, shown by filled 
(empty) symbols. Charge stripes are arranged along the y di- 
rection in all cases. Energies (per site) are reported in units of 
t, for the case t/ J = 3. The right axis is a temperature scale, 
assuming J — 2000/('. The points considered correspond to 
distortions of 10% and 20% in the two directions. Inset: En- 
ergy difference among the two possible stripe arrangements, 
parallel and perpendicular to the tilt axis, for BC and SC 
wavefunctions. Arrows correspond to the predictions given 
by the Hellman-Feynman theorem for the largest system. 



is not accurate in this case. The resulting picture consists 
of two different types of bonds, some associated to the 
hole-rich region, and others where the superconducting 
order parameter, although sUghtly enhanced or depressed 
locally, is very close the one found for the homogeneous 
RVB state. This feature, together with the non-uniform 
charge distribution, has strong similarities with the pic- 
ture which could be derived from the STM experimental 
findings of Refi^. 



where A| = |cj_|_g ^ — j^c^+a t creates a singlet pair 
of electrons among nearest neighbors for each indepen- 
dent site s of the unitary cell in Figure [T] and a is the 
unit vector, that specifies the bond direction (along x or 
y). We extract the pairing order parameter by taking the 
square root Psirmax) of the pairing correlation function 
at the maximum distance r^ax for different cluster 
sizes, for all independent sites. Notice that the VMC 
and RMFT profiles are in agreement. However, within 
VMC, we find, for all different bonds, that the pairing or- 
der parameter is finite, see Figure [3] Remarkably, this is 
the case also for the bond where the variational parame- 
ter Aij = 0, signaling that the Gutzwiller approximation when the Ay = bonds lie along the direction of the 



IV. LATTICE DISTORTION 

Considering the t — J Hamiltonian of Eq. ([T]), lattice 
distortion is set by imposing two different along the x 
and y direction. We find that setting ^ ay further in- 
creases the stabilization of the non-homogeneous states. 
This feature, that was already observed in the context 
of spin-stripesiiii^, extends its validity in the quite dif- 
ferent case of superconducting-stripe states. Considering 
the case of DRVB states, we find that the bond-centered 
stripes gain a remarkable amount of energy upon dis- 
tortion, see Figure [H In particular, we find that the 
superconducting inphase domain stripes are stabilized— 
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FIG. 5: (Color online) Energy difference between the BC- 
ttDRVB (SC-ttDRVB) and the homogeneous RVB stateSI as 
a function of the lattice distortion Ox/ay, shown by filled 
(empty) symbols. Charge stripes are arranged along the y 
direction in all cases. Energies (per site) are reported in units 
of t. The points considered correspond to distortions of 10% 
and 20% in the two directions. 
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FIG. 6: (Color online) t — t' — J model: Variational energy 
(in units of t) vs. t' /t for the uniform RVB (circles), the 
BC-ttDRVB (squares) and the BC-DRVB (triangles) wave- 
functions. 



ROLE OF t' 



tilt axis, where the hopping and exchange parameters are 
the largest. Assuming by convention the hole-stripe fixed 
along the y direction, we indeed found the bond-centered 
hole stripes favored for ax > cty (tilt axis along x, with 
the stripe perpendicular to it) and the site-centered hole 
stripes favored for < ay (tilt axis along y, i.e., par- 
allel to the stripes, similarly to what found for antifer- 
romagnetic stripes in Refiii). The stability of the two 
possible stripe arrangements is evident in the inset of 
Figure m where AE, i.e., the energy of the stripe parallel 
to the tilt axis minus the energy of the stripe perpen- 
dicular to it, shows a different behavior for BC and SC 
stripes. Note that all data are in agreement with cal- 
culations using the Hellman-Feynman theorem (provid- 
ing directly the slopes in the limit a^-^y) — 1), valid for 
small deformations, which gives AE/t = 0.00520 At/t 
for bond-centered stripes and AE/t = — 0.0G167At/t 
for site-centered stripes. Remarkably, the RMFT ap- 
proach also confirms that these superconducting mod- 
ulated states are stabilized by tetragonal distortion (not 
shown). Instead, in the case of 7r-shift domain stripes, 
which involves higher energy scales, there is no appre- 
ciable difference in the behavior of the bond- and site- 
centered geometries with respect to the distortion axis, 
see Figure [5l Indeed, lattice deformation lowers the 
ttDRVB energy in both cases with respect to the uniform 
RVB state, with a slightly more favorable contribution 
when the distortion axis lies on the y direction, i.e. par- 
allel to the stripes. 



Since the energies of the stripe states are very close to 
the homogeneous RVB state, some further contribution 
to the microscopic Hamiltonian could finally stabilize the 
bond-centered hole-stripe states as found in experiments. 
Here, we consider an additional hopping integral among 
next-nearest neighbor copper atoms, which is relevant, 
in many cuprates, for a correct description of their Fermi 
surface. The t — t' — J Hamiltonian is obtained from the 
t — J model of Eq. ([T|) by adding: 



Ht^t'-j^Hu + t' {ctcj, + h.c.) (5) 



where denotes next-nearest neighbors. Here, we 

take i' > and t/J = 3, being relevant for hole-doped 
cuprates. In the presence of i', a further variational 
hopping parameter Xij among next-nearest neighbors is 
added to the mean- field Hamiltonian of Eq. ([2]). Con- 
sidering the uniform RVB state, x'ij is homogeneous for 
all sites. Instead, in the case of stripe states, four differ- 
ent x'ij are optimized, according to the symmetries of the 
ttDRVB and DRVB unitary cells. In Figure [6] we com- 
pare the variational energies of the two candidate stripe 
states, the 7r-phase domain RVB stripe and the inphase 
domain RVB stripe, with the energy of the homogeneous 
RVB state for different t' . It turns out that the DRVB 
stripe is the most stable inhomogeneous state, but, in 
both cases, the effect of t' is irrelevant for the stabiliza- 
tion of superconducting stripes. 
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FIG. 7: (Color online) Energy of the homogeneous RVB state 
as a function of the Coulomb-repulsion strength Vjt^ in the 
presence of a short-range Jastrow factor (stars) and a long 
range Jastrow (squares). Inset: Energy difference between 
the SFP and the homogeneous RVB state as a function of 
the Coulomb-repulsion strength V jt. Both wavefunctions are 
optimized in the presence of a short-range (long-range) Jas- 
trow factor, the energy difference being represented by stars 
(squares). 



VI. ROLE OF NON-LOCAL COULOMB 
REPULSION 

Finally, we consider the effect of a long-range Coulomb 
repulsion, which could play an important role in the STM 
measurements. To this purpose we consider the t — J — V 
model: 



Ht-.]-v 



H. 



(6) 



where n is the average density of electrons. The poten- 
tial has the form of a screened Coulomb repulsion 
_ y exp(-r/io) ^ ^ being the periodized distance be- 
tween different lattice sites. In the following, we have 
fixed la = 4 and considered the properties of the sys- 
tem as a function of V. At the variational level, we find 
that, in the presence of a non-local Coulomb repulsion, 
a long-range Jastrow factor acting on the fully-projected 
mean- field state |^'g} is needed: 



(7) 



where Vj — exp 



2 



with Vi 



ational parameters which depend only on the distance 
\i — j\ and \^g) is defined in Eq. This term takes 
into account the two-body correlations along all possible 
distances, and is necessary to get accurate wavefunctions 
in the presence of K . In Figure [7] we show the en- 
ergy of the homogeneous RVB state, optimized with a 
short-range (i.e., up to nearest-neighbors) and a long- 
range Jastrow factor, respectively. It turns out that, by 



FIG. 8: (Color online) Jastrow parameter vs. distance for 
two tilted lattices on 242 (stars) and 338 sites (circles) and 
different values of V/t. 



increasing V/t, the long-range correlation term plays a 
fundamental role in lowering the energy of the RVB vari- 
ational state. Therefore, the addition of the long-range 
Jastrow term can strongly influence the hierarchy of vari- 
ational phases stabilized at large V/t. Indeed, consider- 
ing relative stability of the SFP with respect to the RVB 
state (see inset of Figure [T]), it turns out that, with a 
short-range Jastrow, the SFP is the lowest-energy state 
already for V ~ 3i, while the addition of a long-range 
Jastrow stabilizes the RVB state up to very large val- 
ues of the Coulomb repulsion. To better understand the 
role played by the Jastrow factor in the presence of V^, 
we optimize the RVB wavefunction for larger clusters. 
In Figure [5] we plot the resulting Jastrow parameters 
in real space. By increasing V , the correlation terms 
Vr increase, as expected. In particular, notice that the 
part at long distances is non- negligible even for r > Iq. 
Since a uniform Jastrow can promote the stability of uni- 
form states and mask possible instabilities towards mod- 
ulated states, we optimize the Jastrow parameters up to 
the next-nearest neighbor distance independently on the 
nonequivalent bonds (i.e. here the depend on i and j). 
The resulting energies are shown in Figure [9] where, even 
though the stripe energies approach the uniform RVB 
wavefunction by increasing V/t, the homogeneous state 
still remains the most stable. Notice that the presence 
of a long-range Jastrow factor stabilizes the BC-ttDRVB 
state w.r.t. the SFP, at V — 0, contrarily to the val- 
ues reported in Table [H However, by increasing V/t the 
SFP wavefunction gains energy and approaches closely 
the uniform RVB state. Moreover, considering the in- 
phase DRVB state, it turns out that for finite V/t the 
energy difference with respect to the RVB state is very 
small, of the order of the error bars. Unfortunately, in 
the presence of a non-local Coulomb repulsion, the main 
limit of this approach corresponds to the necessity of in- 
corporating long-range Jastrow terms in order to reach 
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FIG. 9: (Color online) Energy difference w.r.t. the uniform 
RVB state of the SFP, ttDRVB stripe state (top) and the 
DRVB state (bottom), as a function of V/t. All wavefunctions 
are optimized in the presence of a long-range Jastrow factor. 



a good accuracy. Since the energies into play are very 
small, the restriction of homogeneous Jastrow terms at 
large distances, which governs the low-energy physics, 
could ultimately favor the uniform RVB state and bias 
our variational outcome. In other words, it might not 
be sufficient to assume inhomogeneous Jastrow terms at 
short distances. However, these variational results sug- 



gest already that long-range Coulomb repulsion can be a 
driving force for the stabilization of modulated states. 



VII. SUMMARY AND CONCLUSIONS 

In summary, we have shown that unidirectionally mod- 
ulated superconducting states are remarkably close in en- 
ergy to the uniform RVB state. This is the case both 
for TT-phasc shift domain RVB stripes and for inphase 
domain RVB stripes. Their properties reflect the re- 
cent STM observations, both concerning the modulated 
charge and superconducting features. Besides the lattice 
distortion, long-ranged Coulomb repulsion further stabi- 
lizes the non-uniform superconducting stripe phases. In- 
stead, next-nearest neighbor hopping does not seem to 
play any role in the hierarchy of phases studied here. 
However, although the stripe energies can approach very 
closely that of the homogeneous RVB state, so far the 
uniform solution remains the most stable. In particular, 
we have found that a 7r-phase shift costs more energy 
than the inphase counterpart having simple A = do- 
mains. Experiments suggest that the real stabilization of 
hole-stripes might imply some further factors than those 
investigated here. 

Besides the perturbations to the t — J Hamiltonian 
considered in this paper, it would be interesting to study 
the effect of disorder in HTSC materials which, according 
to recent calculations22,^ could be relevant in stabilizing 
this kind of superstructures. 
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